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ABSTRACT 

With the first phase of the Square Kilometre Array (SKAl) entering into its final pre¬ 
construction phase, we investigate how best to maximise its scientific return. Specif¬ 
ically, we focus on the statistical measurement of the 21 cm power spectrum (PS) 
from the epoch of reionization (EoR) using the low frequency array, SKAl-low. To 
facilitate this investigation we use the recently developed MCMC based EoR analysis 
tool 21CMMC (Greig & Mesinger). In light of the recent 50 per cent cost reduction, 
we consider several different SKA core baseline designs, changing: (i) the number of 
antenna stations; (ii) the number of dipoles per station; and also (iii) the distribu¬ 
tion of baseline lengths. We find that a design with a reduced number of dipoles per 
core station (increased field of view and total number of core stations), together with 
shortened baselines, maximises the recovered EoR signal. With this optimal baseline 
design, we investigate three observing strategies, analysing the trade-off between low¬ 
ering the instrumental thermal noise against increasing the field of view. SKAl-low 
intends to perform a three tiered observing approach, including a deep 100 deg^ at 
1000 h, a medium-deep 1000 deg^ at 100 h and a shallow 10 000 deg^ at 10 h survey. 
We find that the three observing strategies result in comparable (< per cent) con¬ 
straints on our EoR astrophysical parameters. This is contrary to naive predictions 
based purely on the total signal-to-noise, thus highlighting the need to use EoR pa¬ 
rameter constraints as a figure of merit, in order to maximise scientific returns with 
next generation interferometers. 

Key words: galaxies: high-redshift - intergalactic medium - cosmology: theory - 
dark ages, reionization, first stars - diffuse radiation - early Universe 


1 INTRODUCTION 

The gravitational growth and collapse of initial density per¬ 
turbations into the first astrophysical objects (i.e. stars and 
galaxies) signifies the transition from the end of the cosmic 
dark ages and beginning of the epoch of reionization (EoR). 
Emitted ultraviolet radiation from the hrst sources escapes 
into the surrounding intergalactic medium (IGM) forming 
ionized bubbles (H ii regions) in the pervasive neutral hy¬ 
drogen fog which had enshrouded the Universe since recom¬ 
bination. As more astrophysical structures form, these Hii 
regions grow and overlap eventually permeating the entire 
Universe. 

Evidently, the EoR is rich in astrophysical information. 
It provides insight into the formation, growth and evolution 
of structure in the Universe, the nature of the first stars and 


galaxies and their impact on the physics of the IGM, such 
as its ionization state and temperature (see e.g. |Barkana| 


fc Loeb||2007 Loeb fc Furlanetto||2013 Zaroubi||2013 1. The 

most promising cosmological probe of the reionization epoch 
is the radio detection of the redshifted 21 cm spin-flip tran¬ 


sition of neutral hydrogen (see e.g. Furlanetto et al. 2006 
Morales fc Wyithe||2010 Pritchard fc Loeb|2012l. 

In search of the first direct detection from the EoR, sev¬ 
eral dedicated radio interferometers have been constructed 
or are planned. First-generation 21 cm experiments such as 


the Low Frequency Array (LOFAR; van Haarlem et al. 2013 


Yatawatta et al. 2013 

q the Murchison Wide Field Array 

(MWA; Tingay et al.j 

2013f 

and the Precision Array for 


Probing the Epoch of Reionization (PAPER; Parsons et al. 


* E-mail: bradley.greig@sns.it 


^ http://www.lofar.org/ 

^ http://www.mwatelescope.org/ 
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2 B. Greig et al. 


2010[ P]have been operating for the past few years. In princi¬ 
ple, some of these experiments might yield a detection of the 
EoR 21 cm power spectrum (PS) (e.g. a marginal detection 
of the advanced stages of the reionization epoch; [Chapman] 
Mesinger et al.j 20141. Even without a detec- 


et al. 2012 


tion, these experiments enable substantial progress in the 
understanding of various systematics and the characterisa¬ 
tion and mitigation of astrophysical foregrounds crucial for 
future experiments. 

Proposed, second-generatio n experiments such a s the 
Square Kilometre Array (SKA; Mellema et al. 2013|p and 


the Hydrogen E poch of Reionization Array (HERA; Beard- 
sley et al.|2015| p] with substantially increased collecting ar¬ 
eas will almost certainly detect the EoR 21 cm PS. How¬ 
ever, unlike HERA, whose configuration has been adopted 
specifically to maximise its sensitivity to the PS, the SKA 
core baseline design and configuration are yet to be finalised. 
Recently, in the first step to achieving this, the generalised 
design for the SKA was announced. For SKAl-low, the low 
frequency instrument for which the detection of the EoR 
with 21 cm radiation is amongst the primary science cases, 
this allows for 50 per cent of the original outlined desigrj^ an 
approximate total of ~ 130 000 dipole antennae. It is there¬ 
fore crucial to investigate the optimal distribution of these 
dipoles for maximising the available EoR science. 

Several authors have constructed theoretical models to 
predict interferometer sensitivities facilitating the explo¬ 
ration of various aspects of the design optimisation pro- 


cess ( 

Morales & Hewitt 

2004 

McQuinn et al. 

2006] Par- 

sons et al.|2012a||Beards 

ey et al.|20131|Mellema et al.|20131 

Trott 20141. The overall design depends on many factors; 


the size and distribution of the collecting area, the avail¬ 
able frequency coverage, beam size (field of view), tracked 
or drift scanning and balancing the cost of correlating large 


Mellema et al. 

2013 for 

Parsons et al. 

(2012al 


explored the optimisation of an interferometer layout, dis¬ 
cussing the trade-off between improving the raw sensitivity 
through maximising the number of short, redundant base¬ 
lines (reducing the necessary number of correlators) for 21 
cm PS observations relative to the minimum redundant con¬ 
figurations best suited for imaging and foreground charac¬ 
terisation. 

For SKAl-low, we cannot simply assume a maximumly 
redundant design to recover the highest sensitivity measure¬ 
ment of the PS as it is seeking both a statistical character¬ 
isation and tomographic imaging of the reionization epoch 
(Mellema et al. 20131. Furthermore, SKAl-low has addi¬ 
tional science goals beyond the EoB0 including observing 
Hi line absorption, pulsars, magnetised plasma both in the 
galaxy and IGM, radio recombination lines and extrasolar 
planets. Therefore, optimising the final SKAl-low operat¬ 
ing layout requires a more detailed exploration to satisfy all 
science goals, beyond the scope of this work. 


^ http://eor.berkeley.edu/ 
https: //www.skatelescope.org 


® http://reionization.org 
® https://www.skatelescope.org/wp- 

content/uploads/2014/11/SKA-TEL-SKO-0000002-AG-BD- 

DD-Rev01-SKAl_System_Baseline_Design.pdf 

^ http://pos.sissa.it/cgi-bin/reader/conf.cgi?confid=215 


In this work, we focus instead on how best to distribute 
the available dipole antennae specifically targeting the core 
of the assumed layout configuration. Reducing the total 
number of dipoles allows several scenarios to be explored. 
Either (i) the number of dipoles per core antenna station 
can be reduced (increasing the beam size) or (ii) the number 
of core stations can be reduced (retaining the original beam 
size). Additionally (iii) the average baseline length can be 
reduced producing a more compact design, mitigating the 
reduced sensitivity due to the loss of dipoles. However, a 
more compact design may reduce the angular resolution of 
the instrument which could impact on the ability to provide 
images of individual H ii regions in the reionization epoch. 

Another important aspect to consider is the optimal 
observing strategy to maximise the detection of the 21 cm 
PS. The shape and amplitude of the PS varies across the 
reionization epoch and depends on the astrophysics of the 
reionization process. The total noise of a radio interferom¬ 
eter is dominated by the sample (cosmic) variance on large 
scales, while small-scale modes are affected by the intrinsic 
detector (thermal) noise. By increasing the total integration 
time of a single observed patch sky, the thermal noise contri¬ 
bution to the total noise PS can be reduced, increasing the 
sensitivity on small to intermediate scale modes, while the 
contribution from sample variance remains constant. Con¬ 
versely, for N observations of the sky (where the total on 
sky time is fixed and distributed evenly across the N ob¬ 
servations), the contribution from sample variance to the 
total noise PS can be reduced by y/N, thereby increasing 
the sensitivity to large-scale mode^ 

In order to explore the optimisation of SKAl-low, we 
utilise the Monte Carlo Markov Chain (MCMC) based EoR 
analysis tool 21CMMcP] ( jCreig fc Mesinger|2015 1 (hereafter 
GM15). Designed specifically for swift recovery of astrophys¬ 
ical constraints on the EoR given any observed statistical 
measure of the cosmic 21 cm signal, its flexibility is per¬ 
fectly suited for optimising any telescope design and observ¬ 
ing strategy. Importantly, this is the first optimisation study 
to use the precision of recovered EoR model parameters ob¬ 
tained from a mock measurement of the 21 cm signal as 
a figure of merit, improving on previous approaches based 
solely on the total integrated signal to noise (S/N) of the 
measured PS. 

The remainder of this paper is organised as follows. In 
Section]^ we summarise the key ingredients of 21CMMC to 
be used throughout our analysis, before exploring in Sec¬ 
tion]^ how best to distribute the available antenna dipoles 
to maximise the sensitivity of SKAl-low to the PS. With 
our optimal design, in Section]^ we investigate the optimal 
observing strategy which maximises the recovered accuracy 
of EoR astrophysical constraints. Finally, in Section we 
summarise the key conclusions regarding the optimisation of 
SKAl-low and finish with our closing remarks. Throughout 
this work, we adopt the standard set of ACDM cosmologi¬ 
cal parameters: (Hm, Ha, Hb, n, ag, Ho) =(0.27, 0.73, 0.046, 


® In this instance, while the thermal noise per observation will 
have increased, by combining the thermal noise contribution 
across the N observations will provide a slight improvement to 
the small-scale sensitivity. 

® http://homepage.sns.it/mesinger/21CMMC.html 
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0.96, 0.82, 70 kms ^ Mpc ^), measured from the Wilkinson 
Microwave Anisotropy Probe (Bennett et al. 20131 which 


are consistent with the latest results from Planck (Planck 
Collaboration XVI||2014 K 


2 METHODOLOGY 

In GM15, we developed 21CMMC: an MCMC driver 
of the publicly available semi-numerical code 21CMFAST 
( [Mesinger &: Furlanetto| [2007| [Mesinger et al.||2011[ ). We 
showcased how it can be used to recover Bayesian constraints 
on astrophysical EoR parameters, given a mock PS obser¬ 
vation with current and upcoming interferometers. Here we 
focus on the optimisation of SKAl-low, using as a figure 
of merit the precision to which we can recover astrophysical 
EoR parameters from the 21 cm PS. In the remainder of this 
section we summarise the relevant components of 21CMMC 
necessary for performing this work, and defer the reader to 
GM15 for more detailed discussions. 


2.1 Modelling reionization 


Within 21CMMC reionization is modelled using 21CMFAST, 
where 3D realisations of the IGM density, velocity, source 
and ionization helds are generated by evolving and smooth¬ 
ing a cubic volume of the linear density held using the 


Zel’dovich approximation (Zel’dovich 19701. ionization helds 


are then estimated using an excursion-set approach (Furlan- 
|etto et S^|2004[ ), comparing the time-integrated number of 
ionizing photons to the number of baryons within spherical 
regions of decreasing radius, R. Any cell within the simula¬ 
tion volume is classihed as fully ionized if. 


C./'coll (.r, Z, R, Alniin) ^ 1, 


( 1 ) 


where is the ionization efhciency which describes the con¬ 


version of mass into ionizing photons (see Section 2.1.11 
and fcou{x, z, R, Mmin) is the fraction of collapsed matter 
within a spherical radius R residing within haloes larger 
than Mmin (Press fc Schechter|l974 Bond et al.[l991 Lacey 


& Gole 1993| Sheth & Tormen 19991. For cells not fully 
ionized (partial ionization) the ionized fraction is set to 
C/coii(a:, z, Rceii, Mmin), where Rceii is the minimum smooth¬ 


ing scale of the cell. 

In GM15, we showed that 21CMMC can accommodate 
various parametrisations of EoR astrophysics. For simplic¬ 
ity, here we use a popular three parameter model, consist¬ 
ing of: (i) the ionizing efficiency of high -2 galaxies; (ii) the 
mean free path of ionizing photons; (iii) the minimum virial 
temperature hosting star-forming galaxies. Although overly 
simplistic (as galaxy and IGM properties are averaged over 
redshift and/or halo mass dependences), this model provides 
the flexibility to describe a broad range of EoR signals while 
providing a straightforward physical interpretation, which 
we now briefly summarise. 


2.1.1 ionizing efficiency, g 


The ionizing efficiency of high -2 galaxies (Eqn. can be 
expressed as 


c 


30 


/esc \ 

0.3 J 



N-, \ 
4000 J 


1 + Ur. 


( 2 ) 


where, /esc is the fraction of ionizing photons escaping into 
the IGM, /* is the fraction of galactic gas in stars, is the 
number of ionizing photons produced per baryon in stars 
and rirec is the typical number of times a hydrogen atom 
recombines. Although several of these parameters are un¬ 
certain, we showcase some possible values above, resulting 
in our fiducial choice of C) = 30 used for our mock observa¬ 
tion. Here we explore the range G [5,100]. For reference, 
we also show on the top axis of the corresponding panels the 
parameter /esc, converted from / using the fiducial choices 
shown in Eqn. (|^. 


2.1.2 Mean free path of ionizing photons within ionized 
regions, Rmfp 

The propagation of escaping ionizing photons through the 
IGM depends strongly on the abundances and properties of 
the absorption systems (Lyman limit as well as more diffuse 
systems). Below the resolution limits of EoR simulations, 
these act as photon sinks roughly dictating a maximum 
physical scale to which Hii bubbles are capable of growing 
around ionizing galaxies. Typically in EoR modelling this is 
parametrised by a maximum horizon for the ionizing pho¬ 
tons, Rmfp, corresponding to the maximum hltering scale 
in excursion-set EoR models (e.g. [Furlanetto fc Oh||2005p . 
Physically, this is often regarded as the mean free path of 
ionizing photons within ionized regions (e.g. O’Meara et al. 
2007[ [Prochaska, fc Wolf^|2009[ [Songaila fc Cowie||2010| 


McQuinn et al. 20111, though during the EoR this maxi¬ 


mum horizon can be smaller than the instantaneous mean 
free path ( [Sobacchi fc Mesinger|20T^ . Motivated by models 
containing sub-grid physics to model recombinations within 


these small absorption systems (e.g. Sobacchi & Mesinger 
[Mil l, we adopt a fiducial range of Rmfp £ [5, 20] cMpc, and 
the value of Rmfp = 15 Mpc for our mock observation. 


2.1.3 


Minimum virial temperature of star-forming haloes, 

ry-imin 

vir 


We dehne the minimum threshold for a halo hosting a star¬ 
forming galaxy in terms of its virial temperature, TT/“. This 
minimum halo virial temperature is set by the requirement 
of efficient: gas accretion, cooling and retainment of super¬ 
novae outflows. For reference, the virial temperature can be 
related to the halo mass via, (e.g. Barkana fc Loeb|2001 1: 

1 / p /o A \ -1/2 

voe/ 


M 


10*h" 


Tvi 


1.98 X 


IQi K J 


\ 3/2 


V 


IStt^ j 

1-t-z 

10 


- 3/2 


MA¬ 


O’ 


( 3 ) 


where p is the mean molecular weight, = Hm(l + 
2)^/[Oni(l + z)^ -(- Ha], and Ac = 187r^ -|- 82d — 39d^ where 
d = — 1. Here we allow this threshold temperature to 

vary within the range T™" € ]10"^,2 x 10®] K, with the 
lower limit corresponding to the atomic cooling threshold 
and the upper limit corresponding roughly to the observed 


Lyman break galaxy candidates. We choose = 30 000 
K (logT™/" = 4.48) for our mock observation, resulting in 
a reionization history which matches recent estimates based 
on the CMB optical depth (Planck Collaboration XHI 20151. 
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4 B. Greig et al. 


2.2 Telescope noise profiles 

In order to explore the optimal core baseline design and 
observing strategy for SKAl-low, we compute the theoret¬ 
ical noise PS using the publicly available python module 


21CMSENSi ^Pober et al. 


2013 


defer the reader to Parsons et al. (2012a I and Pober et al. 


20141. For further details we 


(2013 20141, and instead briefly summarise the important 


assumptions. 

The thermal noise PS is computed at each Mu-cell ac¬ 
cording to the following (e.g. Morales|2005 McQuinn et al. 
^006}[Pober et al.|2014| ), 


(4) 

where X^Y is a cosmological conversion factor between ob¬ 
serving bandwidth, frequency and comoving volume, fl' is 


a beam-dependent factor derived in Parsons et al. (20141, t 


is the total time spent by all baselines within a particular 
k mode and Pays is the system temperature, the sum of the 
receiver temperature. Tree, and the sky temperature Psky 
Of great concern for 21 cm experiments is the es¬ 
timation and removal of the bright foreground emission. 
Throughout this work, we conservatively choose to adopt 
the ‘moderate’ foreground removal strategy of [Pober et al.| 


(20141. We excise Fourier modes which are localised into 


et af.||2010 

[Morales et af.|2012| Parsons et al.[[2012b| 

Trott 

et al.|2012 

[Vedantham et al.|2012[ Thyagarajan et al. 

20T^ 


foreground ‘wedge’ from the observing window, we fix its 
location to extend Afc|| = 0.1 h Mpc“^ beyond the horizon 
limit ( jPober et al.||2014| |. This foreground removal strategy 
additionally includes the summation over instantaneously 
redundant baselinewhich can be coherently added to re¬ 
duce the uncertainties in the thermal noise estimates ( jPar-j 
sons et al.|2012a|. 


Within 21CMSENSE the total noise PS is then computed 
by performing an inverse-weighted summation over all in¬ 
dividual modes (e.g. jMcQuinn et al.|2006 l, combining both 
the contribution from the sample variance and the thermal 
noise (equation]^. 

- 1/2 


5A 


T+s{k) — [ 


{AlAk) + Xl,{k)y 


( 5 ) 


where {k) is the total uncertainty from thermal noise 

and sample variance in a given fc-mode and A2i(fc) is the 
cosmological 21 cm PS from 21CMFAST. Here, we assume 
Gaussian errors for the cosmic-variance term, which on large 
scales is a good approximation. As in GM15, in addition to 
the aforementioned foreground ‘wedge’, we include another 
level of conservatism to our errors, applying a strict fc-mode 


https://github.com/jpober/21cmSense 

It is important to note that this observing window is still 
not completely free of foregrounds. Fourier transforming the data 
cubes introduces side-lobe noise that can leak into this observing 
window and hence foreground removal is still required. However, 
conveniently the majority of this noise still occurs within this 
foreground ‘wedge’. 

Note, all 21 cm experiments contain redundant baselines, how¬ 
ever, these are not always instantaneously redundant. 


cut, |fc| < 0.15 Mpc“^. Due to the location of the foreground 
‘wedge’ there is limited information available below this k- 
mode, however, modes above can still be affected by the 
foreground ‘wedge’ due to its spherically asymmetric struc¬ 
ture, hence we choose to apply both cuts. 

Previously, both jPober et al.| ( |2014| ) and GM15 assumed 
a drift scan observing mode for the SKA when estimating re¬ 
covered EoR parameter constraints to aid a like-for-like com¬ 
parison between the various 21 cm experiments. However, 
SKAl-low, is designed to operate in track scanning mode. 
In order to accurately discuss the optimisation of SKAl- 
low, we therefore adopt tracked scanning into our analysis. 
The major difference between the two scanning modes is on 
which spatial scales they are able to maximise the sensi¬ 
tivity. Drift scanning experiments are only able to observe 
chosen regions of the sky with low system temperatures for 
short periods of time (the actual time is dependent on the 
beam size) owing to the Earth’s rotation. The short integra¬ 
tion times result in more thermal noise, although they allow 
several independent patches of sky to be observed per ob¬ 
serving run, reducing the sample variance. The opposite is 
true for tracked scanning. As a result, a drift scanning mode 
has improved sensitivity on larger spatial scales (smaller k) 
whereas a tracked scan is favourable on moderate to smaller 
spatial scales. 

Throughout this work we assume for SKAl-low that 
we are able to perform a single 6 h tracked scan on a patch 
of sky per night. This is a conservative choice to ensure our 
tracked regions of sky remain well above the horizon for 
the duration of the integration, mitigating any attenuation 
which may occur from observing close to the horizon. For all 
observations (i.e. our redshift range) we assume a bandwidth 
of 8 MHz and model the frequency dependent sky ternpera- 
ture as Tsky = 60 ( 300 mhz ) ^ ^ (Thompson et al.j2007 1 

while adopting Tree = O.lTsky-l-dO K ior the antenna dipoles. 


2.3 Recovering astrophysical constraints with 
21CMMC 

The large instantaneous bandwidth coverage available with 
the SKA-low (e.g. 50-350 MHz) should facilitate combining 
multiple epoch observations of the 21 cm PS to improve the 
overall constraining power. The improvements arise from the 
effective sampling of the redshift evolution in the 21 cm PS 
on large spatial scales, which are most sensitive to the EoR 
physics (see Fig. 2 of GM15). 

Our mock observation is generated using a 500 Mpe, 
756^ simulation, assuming = 30, i?mfp = 15 cMpc and 
T™" = 3 X 10^ K. All astrophysical constraints are recov¬ 
ered assuming the linear combination of the individual like¬ 
lihoods of the PS observations at 2; = 8, 9 and 10. By assum¬ 
ing an 8 MHz bandwidth, which corresponds to a Xz < 1, 
we ensure there is no overlap in redshift space for these ob¬ 
servations. 
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Design 

layout 

Diameter 

(m) 

Single beam FoV 
(deg^ @ 150 MHz) 

^ Antenna 
stations 

Collecting area 
(m2) 

Description 

#1 

35.0 

14.2 

433 

416 595 

Original SKA design, number of stations halved 

#2 

24.8 

28.4 

866 

416 595 

Original SKA design, dipoles per station halved 

#3 

35.0 

14.2 

433 

416 595 

Baselines reduced by \/2, number of stations halved 

#4 

24.8 

28.4 

866 

416 595 

Baselines reduced by v^, dipoles per station halved 


Table 1. Summary of the different design layouts we consider for SKAl-low. The recently announced general design calls for a 50 per 
cent reduction of the originally planned SKAl-low, therefore we consider various ways of distributing such a reduction. We compute the 

single beam field of view (FoV) using FoVbeam = j —) ’ 


3 OPTIMAL CORE BASELINE DESIGN 

The original design brief for SKAl-low core consisted of 
~250 000 dipole antenna distributed across 866 cor an¬ 
tenna stations. Each station is 35 metres in diameter and 
arranged into a densely configured compact core contain¬ 
ing ~75 per cent of the antenna stations out to ~1 km. 
Depending on the observing frequency each station would 
have a field of view of ~20 deg^. With a 50 per cent reduc¬ 
tion to the original SKAl-low design, the question is, how 
best can these antenna dipoles be distributed to maximise 
the scientific returns for the EoR? 

There are several options which can be considered. To 
reduce cost one can either: (i) cut the total number of core 
antenna stations; or (ii) cut the number of dipoles per sta¬ 
tion (subsequently increasing the field of view). Moreover 
for both (i) and (ii) one could also reduce the average base¬ 
line lengths for the core antenna stationj^ increasing the 
sensitivity over a narrower fe-range. The reduction in the 
spatial resolution which comes with the compactification of 
the core, does not impact our EoR parameter constraints, as 
the bulk of the constraining power comes from larger scales. 
Nevertheless, the loss in spatial resolution would adversely 
impact the proposed direct imaging of the EoR. It is still not 
quantifiably clear what we will learn from such images. A 
more detailed exploration, beyond the scope of this present 
work, is required to quantify how an image-based EoR anal¬ 
ysis is impacted by the choice of baseline lengths. 


3.1 SKA layouts 

Prior to exploring these options, we must first describe the 
layout with which we choose to model SKAl-low. It is 
planned that SKAl-low will approximate a circularly sym¬ 
metric array whose stations are randomly distributed and 
fall off with radius approximating a Gaussian distribution. 
At the same time, it should be as compact as possible, with 


Here, we refer only to the core of SKAl-low. A further 45 
stations were to be configured into spiral arms extending outward 
as far as 50 km, to allow maximum baseline lengths of ~100 km. 
We do not consider these stations as they add very little to the 
sensitivity of the 21 cm PS. Instead these aid calibration and 
foreground removal ||Mellema et al.|2013||. 

Note, we consider only reducing the average baseline lengths 
of the core antenna stations, baselines beyond ~10 km used for 
calibration and foregrounds would remain in their original con¬ 
figuration. 


the central core requiring a tightly packed configuration with 
a filling factor as close to unity as possible. The tightest 
packed configuration for circular stations is the hexagonal 
configuration. Indeed, Parsons et al. (2012a I showed that a 
hexagonal configuration provides the greatest level of max¬ 
imum baseline redundancy increasing the overall telescope 
sensitivity to these Fourier modes, which in turn are best 
suited for the statistical recovery of the 21 cm PS. 

Although this configuration is by no means the final 
baseline design, for the remainder of this work, we as¬ 
sume that SKAl-low has a central, densely packed core 
in a hexagonal configuration out to a certain radius, be¬ 
yond which the remaining stations are randomly distributed 
falling off with radius (neglecting stations in the spiral arms 
for calibration and foreground removal). The radius of the 
hexagonal core is not fixed, but chosen per design to ap¬ 
proximate as close as possible the required distribution of 
stations with radius and maintaining a filling factor as close 
to unity as possible, in general accordance with the SKA 
baseline design. We shall consider four distinct designs for 
the distribution of the dipole antennas, described below and 
summarised in Table [T] 


• Design ^1: Original design for SKAl-low, halving 
the number of antenna stations down to a total of 433 
(retaining the original number of dipoles per station, i.e. 
the original field of view). This design retains 75 per cent 
of the stations within a radius of ~1 km, and contains 61 
stations in a hexagonal configuration in the inner core. 

• Design 7)^2: Original design for SKAl-low, halving the 
number of dipole antennas per station. We retain the origi¬ 
nally planned 866 stations, with a reduction in the station 
diameter of y/2 (doubling the field of view of the antenna 
stations). To construct this design, we originally assume 
an antenna station diameter of 35 m before shrinking the 
station diameter. A total of 127 stations are included in the 
hexagonal inner core, and 75 per cent of the stations are 
within a radius of ~1 km. 

• Design ^3: Same as design ^1, except we additionally 
reduce the average baseline length by \/2 producing a 
more compact overall design. With this, we increase the 
size of the hexagonal inner core to 91 stations, and now 
contain 75 per cent of the stations within a radius of ~700 m. 

• Design 7)^4: Same as Design ^2, reducing the average 
baseline length by \/2 to produce a more compact overall 
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Figure 1. A schematic of the four core baseline designs we consider for SKAl—low. Design ^1 (top left): 433 core antenna stations 
with a diameter of 35 m. Design ^^2 (top right): 866 core antenna stations with half the number of dipoles per station as design 
(station diameter 35/\/2 m). Design ^3 (bottom left): Same as design except reducing the baseline length by y/2. Design (bottom 
right): Same as design ^^2 except reducing the baseline length by \/2. The relative physical size of the filled circles correspond to the 
true diameter of the antenna stations. For all designs the core antenna stations are normally distributed with radius, including a central 
densely packed ‘inner’ core in a hexagonal configuration (red circles). For design ^1 and ^^2, blue, green and yellow circles denote stations 
within a 600 m radii (blue dashed curve), 1000 m radii (green dashed curve) and stations beyond 1000 m from the core respectively. For 
design ^3 and these radii are reduced by \/2- 


design. With the increased number of stations in this design 
(866), the inner hexagonal core retains the number of sta¬ 
tions as design ^^2 but is shrunk to be fully compact. This 
design maintains 75 per cent of the stations within a radius 
of -700 m. 

For a visual description of these designs we provide 
Fig-H All filled circles correspond to the relative diameter of 
the antenna stations. Within the inner hexagonal core sta¬ 
tions are denoted by red circles, while blue, green and yellow 


circles correspond to stations within certain radii (empha¬ 
sised by the circular dashed lines). For designs and #2, 
these correspond to stations within a 600 m radii, 1 km radii 
and stations outside a 1 km radii of the centre, respectively. 
For designs ^3 and #4 (with reduced baselines), these radii 
are reduced by \/2. In Fig. we additionally provide the 
cumulative number of antenna stations with radius to bet¬ 
ter highlight the distribution of antenna stations amongst 
our four chosen designs. 
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Figure 2. The cumulative distribution of SKAl-low telescope 
antenna as a function of radius for all instrument designs we con¬ 
sider. The red curves correspond to retaining the original number 
of core antenna stations by halving the number of dipoles per core 
antenna station, while the yellow curves correspond to halving the 
number of core stations. Dashed curves correspond to the original 
baseline lengths while the solid curves denote designs where the 
average baseline length is reduced by \/2- 


3.2 Power spectrum sensitivity per beam 

In Fig. we compute the total noise PS (Eqn. from 
21CMSENSE for each design assuming a single beam point¬ 
ing integrated for a total of 1000 h. Here, and throughout 
the remainder of this work, we assume that SKAl-low can 
always be calibrated to reach its theoretical noise limit as 
defined by the total noise PS. We show the total theoreti¬ 
cal noise PS for all designs at each redshift for our multiple 
epoch observing strategy (a = 8, 9 and 10, left, centre and 
right panel, respectively), where the solid black curve cor¬ 
responds to the mock 21 cm PS. The red and yellow curves 
denote halving the number of dipoles per station or halv¬ 
ing the total number of stations, respectively. The dashed 
curves denote the retention of the original baseline distribu¬ 
tion, whereas the solid curves denote reducing the average 
baseline lengths. Dot-dashed and dotted curves denote the 
contribution from only the thermal noise component for the 
dashed and solid curves respectively. 

The major difference between the curves in Fig. [^sterns 
from the increased held of view of designs #2 and #4, rel¬ 
ative to the other designs. Reducing the number of dipoles 
per core antenna station, doubles the effective held of view 
(see Table [^, improving the sensitivity to Fourier modes 
most relevant for measuring the 21 cm PS. The larger ehec- 
tive sky area observed per beam pointing reduces the error 
contribution from sample variance, as well as that of the 
thermal noise (see Eqn. since each mode in the sky can 
be observed with the same S/N per resolution element over 
a larger part of the sky. As a result, for the same total ob¬ 
serving time designs and #4 are optimal. 


3.3 Power spectrum sensitivity per fixed survey 
area 

In order to highlight this, in Fig. we provide the same 
as Fig. except we normalise the sensitivity curves to the 
same total survey area footprint. Here, we consider a total 
area of 100 deg^ on the sky observed at 1000 h per beam. 
Note for this survey, designs ^1 and ^3, with half the held 
of view, ehectively take twice as long in terms of total time 
spent observing on the sky, due to the increased pointings 
required to obtain the hxed survey area. With the impact of 
held of view now normalised out, we hud designs #1 and #2 
to be equivalent, and the same for designs #3 and ^4. This 
is theoretically expected as by observing the same total sur¬ 
vey area, all designs have the same error contribution from 
sample variance. Additionally, Fig. [^clearly emphasises the 
impact of reducing the average baseline lengths to produce 
a more compact instrument. On large scales there is no im¬ 
provement, as the total error is completely dominated by 
sample variance, however, on intermediate to smaller scales 
the impact becomes more pronounced. By reducing the aver¬ 
age baseline length, in effect we increase the number of core 
antenna stations separated by baselines sensitive to these 
Fourier modes, improving the overall sensitivity of the in¬ 
strument. However, if we were to continue down to even 
smaller scales we would eventually see a turn over, with de¬ 
creased sensitivity for very small scale modes. 

Clearly, from Figs|^and|^the optimal design for SKAl- 
low is then design ^4, from which we beneht both from the 
reduction in baseline lengths, and the decreased total on sky 
observing time (due to a larger FoV). However, our choice 
of design ^4 over design ^(^3 does not factor in other con¬ 
siderations. For example, with double the number of core 
antenna stations than design ^3, design ^4 will have a sub¬ 
stantially higher computational cost for correlating the sig¬ 
nal. Additionally, construction costs will be larger owing to 
more stations being required. Despite these caveats, for the 
remainder of this work, we shall consider design ^4 to be 
the optimal design for the SKAl-low. Moreover, we will only 
consider surveys of total fixed area and not total observing 
time. Therefore, any conclusions drawn would be equivalent 
between designs #3 and #4. However, if for example the 
observing time was additionally fixed (owing to survey re¬ 
strictions) , design #4 would outperform design #3 by effec¬ 
tively being able to spend twice as long per beam pointing, 
providing a gain in S/N over the entire survey area (\/2 for 
imaging and factor of two for the PS sensitivity). 


4 OPTIMAL OBSERVING STRATEGY AND 
EOR PARAMETER FORECASTS 

With our optimal design for SKAl-low, we now seek to de¬ 
termine the optimal observing strategy. The overall sensitiv¬ 
ity of a radio interferometer to the 21 cm PS is dictated by 
the relative contributions from the sample (cosmic) variance 
and thermal noise. On large scales, the noise PS is dominated 
by the contribution from sample variance which scales as 
Agy oc an assumed fixed total collecting area of 

the instrument; e.g. Mellema et al.|[2013 ), where Ha is the 
total survey area. On small scales, the noise PS is instead 
dominated by thermal noise and scales as Air OC Da^^^A> 
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Figure 3. The estimated total noise PS for a single beam pointing of 1000 h for each baseline design we consider for SKAl-low at three 
different observational epochs (2 = 8, 9 and 10). The solid black curve corresponds to our fiducial reionization model, with (i^, Rmfp^ 

= (30, 15 Mpc, 3 X 10*^ K) with a corresponding neutral fraction of nni = 0.5, 0.7 and 0.8, respectively. Red and yellow curves denote 
designs whereby we either halve the number of dipoles per core antenna station or the total number of antenna stations, respectively. 
Dashed curves denote the retention of the original baseline lengths, whereas the solid curves correspond to reducing the average baseline 
length by y/2 producing a more compact instrument. Dot-dashed and dotted curves indicate the thermal noise contribution only for the 
dashed and solid curves, respectively. Vertical dashed lines at A: = 0.15 and 1.0 Mpc“^ demarcate the physical scales within which we fit 
the 21cm PS within 21CMMC. 
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Figure 4. Same as Fig. except now we compare our four baseline designs for the SKAl—low for a survey with a total observed area 
of sky of 100 deg^ at 1000 h. This, as expected, mitigates the differences between the two options of halving the number of stations or 
halving the dipoles per station and emphasises the improvement gained on moderate to small scales by reducing the baseline lengths. 
Note, the total time on sky for designs #1 and #3 is twice as long as for designs #2 and #4 due to the reduced field of view (see Table[^. 


where t is the time spent on each pointing. To maximise the 
EoR science return, would a fixed amount of telescope time 
be better spent increasing Q,a or fiR] 

SKAl-low anticipates performing a three tiered obser¬ 
vational strategy in order to balance between the compet¬ 
ing noise contributions, resulting in a relatively uniform S/N 
coverage over a broad stretch of fc-modes. These inclnde a 


deep, medium-deep and a shallow survey (Koopmans et al. 
20151 which we outline below; 


of a 100 deg patch of sky observed at 1000 h. Such a 
survey is designed for the imaging of the reionization epoch, 
seeking the characterisation and detection of individual 
Hii regions. This approach minimises the thermal noise 
contribution to the total noise PS, favouring intermediate to 
small-scales of the 21 cm PS at the expense of an increasing 
sample variance component reducing the sensitivity to the 
large-scale modes in the PS from this survey (note however 
that the thermal noise as shown in Fig. [^and[^is smallest 
on the largest angular scales; see e.g. Koopmans et ah] 


2 '-'11 LllC ICUlgjCOL CLllgLllCtl D^.^<X1C3, OCC C.g. 

100 deg at 1000 h (deep): a single, deep integration 20151. This survey would also extend furthest in redshift, 


allowing us to probe the pre-reionization epoch. 


In obtaining these scaling relations, we have adopted our op¬ 
timal SKA design, design #4, which has a fixed field of view. 
Therefore to observe the requisite total sky area we assume it 
can always be obtained by tiling the survey by N (~ flA/f^Fov) 
observations of the antenna station field of view for t hours. 


• 1000 deg^ at 100 h (medium-deep): an intermediate 
observing strategy, balancing sample variance and thermal 
noise. 
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• 10 000 deg^ at 10 h (shallow): a short integration time 
per pointing^to build a wide/shallow survey at 10 000 deg^ 
after ~ 400 pointings. This survey substantially reduces 
the sample variance contribution, focussing the sensitivity 
of SKAl-low to the largest-scale modes at the expense of 
an increasing thermal noise contribution on all other scales. 

It is important to note that we do not seek to choose 
only one of the above observing strategies as SKAl-low aims 
to perform all three in order to obtain a uniform PS S/N cov¬ 
erage over the widest possible fc-range. Instead, we provide 
forecasts for how well each observing strategy will recover 
the astrophysical parameters of the EoR. 

In Fig. we present the total noise PS for each of the 
three observing strategies we consider for the same three 
redshift observations (« = 8, 9 and 10). Again, the solid 
black curve denotes the theoretical 21 cm PS for our fidu¬ 
cial EoR model. The blue, red and yellow curves correspond 
to the three observing strategies, 100 deg^ at 1000 h, 1000 
deg^ at 100 h and 10 000 deg^ at 10 h, respectively, while 
the dashed curves correspond to the thermal noise contribu¬ 
tion. This figure clearly shows how each observing strategy 
benefits from their respective choice of minimising the ther¬ 
mal noise or sample variance contribution. Across all three 
panels, we observe a shift in the minima of the total noise 
PS to larger spatial scales (small k) as the sample variance 
contribution is decreased by shortening the observing time 
per patch of skyp^ The converse is true on small scales due 
to the increasing thermal noise contribution. Together, the 
three observing strategies provide relatively uniform S/N 
coverage over this range of k. 

The total integrated S/N of the 21 cm PS (which is 
EoR model dependent), summed across all Fourier modes, 
is relatively similar across the three observing strategies, 
with the 1000 h observing strategy yielding the highest total 
sensitivity. However, the differences between EoR models are 
greatest on the largest scales shown in Fig. 2 of GM15; we 
therefore expect that the integrated S/N is not a good proxy 
for how well EoR astrophysics can be constrained. 

4.1 Theoretical limit on recovered EoR 
parameters 

As discussed in GM15, even if our EoR parametrisation is 
correct, modelling errors are introduced when simulating the 
EoR with 21CMFAST (compared with an idealised, infinite 
resolution, first-principle numerical simulation). Character¬ 
isation of these uncertainties over a broad parameter range 
(e.g. by calibrating to more detailed radiative transfer simu¬ 
lations) will be essential to maximising the achievable scien¬ 
tific gains from the SKA. Here we begin by forecasting the 

Note, 10 hrs with SKAl-low is roughly equivalent to 600 hrs 
with LOFAR, which aims for a statistical detection of the EoR 
signal l |Harker et al.|[2010| |. Therefore, a single 10 hr pointing 
with SKAl-low will be equivalent to the deepest integrations with 
LOFAR to date. 

Note that at 2 : = 10, and to a lesser extent at z = 9, the total 
noise PS for our shallow observing strategy is not able to reach its 
minima due to the assumed contamination by foregrounds. As a 
result, this will impact on the total available constraining power 
from this strategy as we discuss below. 


recovered EoR parameters without including any additional 
modelling uncertainty, before relaxing this assumption in 
Section 14.21 

In Eig. we provide both the 1 and 2D marginalised 
probability distribution functions (PDEs) recovered from 
21CMMC for our EoR model parameters for the three differ¬ 
ent observing strategies. In this figure, and for the remainder 
of this work we assume uniform (flat) priors on all recovered 
EoR parameters within their allowed ranges outlined in Sec¬ 
tion |2.1[ Diagonal panels show the recovered marginalised 
ID PDFs for each corresponding EoR parameter with the 
vertical dashed line indicating the location of our fiducial 
input value. In the lower left corner panels, we provide the 
1 and 2(7 2D marginalised likelihood contours for our EoR 
model parameters, where crosses denote their fiducial values. 

The confidence limits from this figure are quantitatively 
summarised in the top half of Table We report the me¬ 
dian recovered values, and the 16th and 84th percentiles. 
Assuming the recovered ID marginalised PDFs for each 
EoR parameter can be modelled by a Gaussian distribu- 
tio 43 we find the fractional uncertainty to which we re¬ 
cover these to be 3.7/3.4/3.6 per cent for 9.6/4.1/2.7 per 
cent for Rmfp and 0.5/0.5/0.6 per cent for logio(T“r‘“) for 
our deep, medium-deep and shallow observing strategies, re¬ 
spectively]^ 

All three surveys recover tight constraints on our EoR 
parameters, highlighting the strength of SKAl-low at de¬ 
tecting the 21 cm PS. However, the medium-deep and shal¬ 
low surveys somewhat outperform the deep observing strat¬ 
egy. This, however, is not surprising given that the majority 
of the constraining power on the astrophysical EoR param¬ 
eters arises from sampling the large-scale modes at multiple 
redshifts (GM15). This is especially true for the mean free 
path parameter, which at Rmfp > 15 Mpc does not strongly 
impact the reionization history, but does change the 21 cm 
large scale power (see Fig. 2 in GM15); this signature is 
missed by the deep survey, and thus the resulting contours 
on Rmfp are broader than those corresponding to the other 
surveys. 

It is additionally worth noting that we are only using 
the PS to compute the likelihood. The 21 cm signal during 
the EoR is highly non-Gaussian. Thus higher order statis¬ 
tics, taking advantage of the increased small-scale sensitiv¬ 
ity of the deep survey, could provide further information. In 
principle, the 21CMMC framework allows for arbitrary fig¬ 
ures of merit; we defer such studies of alternative statistics 
to future work. 

Finally, it is important to note that given all three ob¬ 
serving strategies are expected to be performed by SKAl- 
low, we can combine our constraints across all three surveys 
into one single measurement. In doing so, we find improved 

By performing a MCMC recovery of the astrophysical EoR 
parameters, the recovered marginalised PDEs for each parameter 
need not be symmetric. However, to a good approximation we find 
the PDFs of our EoR parameters to closely resemble a Gaussian 
distribution. 

It would be instructive to compare our recovered constraints 
from our this updated SKAl-low compared to the original, larger 
design considered in GM15. However, it is not straightforward to 
do so, as previously we only assumed a single beam pointing at 
1000 h, relative to the fixed survey areas we consider in this work. 
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Figure 5. The estimated total noise PS for the three EoR observing strategies to be performed by SKAl-low. A single, deep 1000 hr 
integration over a 100 deg^ survey footprint (blue curve), an intermediate survey area of 1000 deg^ at 100 h (red curve) and many, short 
10 h integrations to produce a wide/shallow 10 000 deg^ (yellow curve). Dashed curves correspond to the contribution from the thermal 
noise PS only. While shortening the integration time of an individual patch of the sky increases the thermal noise resulting in a reduction 
of the sensitivity on intermediate to small scales, it reduces the sample variance contribution to the total noise improving the sensitivity 
on larger physical scales. 


Observing strategy 


Parameter 



a^Hi 


(with/without modelling uncertainty) 

c 

-^mfp 

logic (T“-) 

2 = 8 

2 = 9 

2 = 10 

100 deg2 @ 1000 h (without) 

30.66/1;?? 

15.94/1;!? 

^•^^-0.02 

0 4S+0-01 

0 71 +0 01 
'-'•'-‘--0.01 

0.84/0;0i 

1000 deg^ @ 100 h (without) 

30.25/1;°? 

15.43/°;gi 

4 4«+0-02 
^•^^-0.02 


'-'•''^-0.01 

0.83/0;0i 

10 000 deg^ @ 10 h (without) 

28.71/0;°° 

14.22/0;?! 

4 4*3+0.02 
^•^'^-0.02 

0 47+0-01 

U.Dy_Q 0;L 

O.OZ_Q Q;L 

100 deg^ @ 1000 h (with 10 per cent) 

30.68/i;i| 

15.49/?;!! 

4 4Q+0-05 

4-4y_o.05 

n 4S“t0.02 

b-4o_o.o2 

n 71 +0 01 

'-'•'-‘--0.01 

n S4+0 01 

U-04_Q Q;L 

1000 deg^ @ 100 h (with 10 per cent) 

30.62/2;°8 

15.12/1;°! 

4 4Q+0-06 
^•^^-0.06 

'^•^^-0.02 

0 71 +0 01 
'-'•'-‘--0.01 

n S4+0 01 

U.04_Q Q;L 

10 000 deg^ @ 10 h (with 10 per cent) 

so.Yomt 

14.96/?;0| 

4 4Q+0-07 

^•^^-0.07 

n 4S“t0.02 

b-4o_o.o2 

D 71 +0 02 

u. /1_0.02 

n S4+0 01 

U.04_Q Q;L 

100 deg^ @ 1000 h (with 25 per cent) 

31.68/|;2^ 

14.8l/i;°o 

4 'SI +011 
^•^-*--0.11 

^•^^-0.04 

D 71 +0 03 

U. /1_0.03 

n R/i+0 02 

b-o4_Q Q2 

1000 deg^ @ 100 h (with 25 per cent) 

31.84/|;00 

14.87/!;°° 

4 “SI +011 
^•^-*--0.11 

^•^^-0.04 

D 71 +0 03 

U. /1_0.03 

n R/1+0 02 

b-o4_Q Q2 

10 000 deg^ @ 10 h (with 25 per cent) 

32.10/S;®? 

14.81/!;°? 

^•^^-0.12 

n 49+^-®^ 
^•^^-0.04 

0 71 

'-'•'-‘--0.03 

n R/1+0 02 
^•°^-0.02 


Table 2. Summary of the median recovered values (and associated 16th and 84th percentile errors) for our three EoR model parameters, 
C: ^mfp S'lid associated IGM neutral fraction, xjji for all considered observing strategies (with a 10 and 25 per cent 

modelling uncertainty and without a modelling uncertainty). Our fiducial parameter set is (^, Rnifpi logio'^yp") = (30, 15 Mpc, 4.48) 
which results in an IGM neutral fraction of xhi = 0.48, 0.71, 0.83 at 2 = 8, 9 and 10 respectively. 


fractional uncertainties on our recover EoR parameters of 
2.1 per cent for 2.2 per cent for i?mfp and 0.3 per cent for 
lotr 


4.2 Recovered EoR parameters including 
modelling nncertainty 

Our limits in the previous section assumed that we are able 
to perfectly simulate the EoR without any numerical errors. 
However, this is certainly not the case. Even large, numeri¬ 
cal simulations differing only in their implementation of ra¬ 
diative transfer can result in notable differences for a fixed 
source model (e.g. Iliev et ^|2006 |. Furthermore, the sim¬ 
plifications and approximations within semi-numerical ap¬ 
proaches result in even larger modelling uncertainties, which 
appear to be of the order of lO’s of per cent on scales rele¬ 
vant for recovery of EoR astrophysical parameters from the 
PS (e.g. Zahn et al.|2011 |. 


Therefore, in this section we present parameter forecasts 
in the presence of modelling uncertainties. Note, however, 
that this uncertainty only arises due to existing differences 
between current simulations and theory, and can therefore 
be reduced and/or statistically characterised in future work. 
For simplicity here we include this modelling uncertainty 
simply as a constant multiplicative factor of the modelled 
21 cm PS, summing it in quadrature with the total noise 
PS (equation]^. Although in principle this should be a cor¬ 
related error, this approach serves to showcase the relevant 
trends, before a more comprehensive statistical characteri¬ 
sation becomes available. 

In Fig. we provide the same 1 and 2D marginalised 
likelihood contours from 21GMMC for our EoR model pa¬ 
rameters, but with the additional inclusion of a 25 per cent 
modelling uncertainty. We see a notable increase in all con¬ 
tour widths. Moreover, all three strategies perform compa¬ 
rably. This highlights the need for parameter studies char¬ 
acterising modelling uncertainty, which is the limiting un- 
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Figure 6. The recovered constraints on our reionization model parameters for three different observing strategies, combining three 
independent (z = 8,9 and 10) observations of the 21 cm PS. We consider a single, deep 100 deg^ observation at 1000 h (blue), an 
intermediate medium-deep 1000 deg^ observation at 100 h (red) and a wide/shallow observation of 10 000 deg^ at 10 h (yellow). Across 
the diagonal panels we provide the ID marginalised PDFs for each of our reionization parameters [f, iJmfp and logic(T™,!") respectively] 
and highlight our fiducial model parameter by a vertical dashed line. Additionally, we provide the 1 (thick) and 2cr (thin) 2D joint 
marginalised likelihood contours for our three reionization parameters in the lower left corner panels of the figure (crosses mark our 
fiducial reionization parameters). 


certainty for all of the SKAl-low surveys. We again caution 
that this only applies to the PS of the signal: additional 
information from non-Gaussian statistics (e.g. a simultane¬ 
ous imaging analysis with the deep survey) could further 
improve on these constraints. 

In the lower half of Tablej^we quantitatively summarise 
the above behaviour. The fractional precision to which we 
are able to recover our EoR parameters for our 25 (10) per 
cent modelling uncertainty is 17.9/18.0/20.2 (7.6/8.2/10.4) 
per cent for 18.5/18.2/18.4 (12.6/11.6/12.0) per cent for 
Rmfp and 2.4/2.4/2.7 (1.1/1.2/1.6) per cent for logigT™/" for 
our deep, medium-deep and shallow observing strategies, re¬ 
spectively. We find both the deep (100 deg^ at 1000 h) and 
medium-deep (1000 deg^ at 100 h) surveys to be essentially 
equivalent for both our adopted choices of modelling uncer¬ 
tainty, with the shallow survey (1000 deg^ at 100 h) per¬ 
forming marginally worse (as noted by an increase in errors 
for both g and logj^gT™/"). 

The source of this shift in preference to deeper surveys 
is straightforward to understand. The fractional precision of 


the astrophysical constraints are strongly dependent on the 
sensitivity to which the large-scale modes are able to be sam¬ 
pled. By including a constant multiplicative error, we now 
substantially increase the errors on large scales, completely 
overwhelming any potential gains by shortening the observ¬ 
ing times (decreasing the cosmic variance). In effect, all ob¬ 
serving strategies now have the same sensitivity to the 21 cm 
PS on large scales. On smaller scales, the smaller thermal 
noise contribution from the deep and medium-deep surveys 
gives them the edge in parameter recovery over the shallow 
survey. However, as noted in GM15, increased sensitivity on 
small scales does not significantly aid EoR constraints across 
multiple epoch observations, as the reionization history is 
still adequately sampled from the large-scales. 

Finally, as in Section [4.1[ with SKAl-low we can com¬ 
bine all three observing strategies to provide improved over¬ 
all constraints on our EoR parameters. In the case of our 
25 (10) per cent modelling uncertainty, we then find the im¬ 
proved constraints of 10.5 (4.9) per cent on 10.7 (6.9) per 
cent on Rmfp and 1.4 (0.7) per cent on logjoF/u"- 


© 0000 RAS, MNRAS 000, 000-000 
































12 


B. Greig et al. 


fesc 



C Rmfj, (Mpc) [K]) 


Figure 7. Same as Fig. [^except now we include a 25 per cent modelling uncertainty of the sampled 21 cm PS within 21CMMC. 
This accounts for the expected fractional variation between the various numerical approximations adopted within the semi-numerical 
simulation code 21CMFAST relative to fully numerical reionization simulations. 


5 CONCLUSION 


The reionization epoch is rich in astrophysics: probing the 
formation and evolution of the first galaxies and their im¬ 
pact on the IGM. Upcoming radio interferometers, such as 
SKA 1-low and HERA, will provide a clear window into these 
processes, using the redshifted 21 cm line. It is therefore im¬ 
portant to optimise their designs and observing strategies, 
so that we can maximise the corresponding EoR science re¬ 
turns. 


Using the MCMC based EoR analysis tool 21CMMC 
I Greig fc Mesinger|2015 1, here we explored the optimisation 
of SKAl-low. Given the recent 50 per cent cost reduction, 
we explored how best to distribute the available resources, 
considering: (i) halving the original number of core antenna 
stations (ii) halving the original number of dipoles per core 
station (maintaining the original number of antenna stations 
and doubling the held of view) and (iii) reducing the average 
baseline length by %/2- We found the optimal core baseline 
design for SKAl-low to be one where the number of dipoles 
per station is halved with a concurrent reduction in baseline 
length. The former reduces the contribution from cosmic 
(sample) variance on large scales by observing a larger area 
of sky for a hxed survey duration, while the latter reduces 


the thermal noise on the relevant fc-modes, resulting in the 
optimal sensitivity across all scales. 


We then explored which planned observing strategy for 
SKAl-low provides the strongest constraints on our EoR 
model parameters from the 21 cm PS. SKAl-low plans to 


perfo rm a three tiered observing approach | Koopmans et al. 
20151 (i) a single, deep 100 deg^ survey area observed at 


1000 h, (ii) a medium-deep survey of 1000 deg^ at 100 h and 
(iii) a wide/shallow survey observing 10 000 deg^ at 10 h. We 
found the shallow survey to recover the tightest constraints 
on our model EoR parameters. By shortening the observa¬ 
tion per field of view (large number of independent helds), 
the sample variance is considerably reduced, increasing the 
overall sensitivity on the largest scales which exhibit the 
largest differences between EoR models. It is important to 
note that this conclusion is contrary to what would naively 
be expected from a simple consideration of the total S/N, 
which is maximised in the deep survey. This highlights the 
need to use EoR parameter constraints as a hgure of merit, 
instead of just the total integrated S/N. 


However, once we included modelling uncertainty into 
our analysis, to account for approximations and simplifica¬ 
tions inherent to our simulations, we found the strategies 
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to be comparable, with the deep survey preforming slightly 
better than the others. Modelling uncertainty dominates on 
large scales for all surveys, thus further motivating param¬ 
eter studies which statistically characterise modelling error, 
in order to maximise the science returns from second gener¬ 
ation interferometers. 

It is worth emphasising that this optimisation study was 
based solely the 21 cm PS as a likelihood estimator. While 
this is the primary EoR science data product of SKAl-low, 
other alternative statistics of the 21 cm EoR signal could 
yield additional information, especially in the deep survey 
where 21 cm structures can be directly imaged. These alter¬ 
native statistics, to be explored in future work, will likely 
benefit from the increased sensitivities on smaller spacial 
scales. 

Finally, we stress that in this study we only focused on 
the EoR. The pre-reionization epoch of IGM heating encodes 
an additional wealth of information about high-energy pro¬ 
cesses inside the first galaxies. For such Cosmic Dawn stud¬ 
ies, deep integrations are needed to push out to the highest 
redshifts possible. 
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